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ps) ■ Abstract 

C ' We prove by explicit calculation that Feynman graphs in noncommutative Yang- 

1-^ ■ Mills theory made of repeated insertions into itself of arbitrarily many one-loop ghost 

, propagator corrections are renormalizable by local counterterms. This provides a 

CN I strong support for the renormalizability conjecture of that model. 

T— I ! 

■ 1 Introduction 

(N ; 

^ ■ It is now commonly admitted that our current concepts about space and time have to be 
y—i . changed when exploring space-time at a very small scale. Indeed, one can show that it is 
O I impossible to locate a particle with an arbitrarily small uncertainty when taking both into 
Q I account the principles of quantum mechanics and general relativity Roughly speaking, 
^ ' one can say that measurements of coordinates on space-time are subject to uncertainty 
• relations, thus ruining all geometrical concepts that have proved to be a guidance principle 
Oh. in elaborating many physical theories. 

r-| ! Following the previously alluded analogy with quantum mechanics, one can try to solve 

^ I this puzzle by assuming that the coordinates themselves are noncommuting objects. Thus, 
the natural extension of geometrical ideas to this new type of coordinates has been called 
"noncommutative geometry". Following even closer the ideas and methods of quantum 
mechanics, we are led to assume that these noncommutative coordinates are represented as 
a subalgebra of the algebra of operators acting on a Hilbert space. This is the framework 
of the theory pioneered by A. Connes 0, which allows us to make use of the powerful 
tools of functional analysis. Within this framework, an analogue of gauge theory has been 
developed, even with non trivial topological properties, and it has already proved to be 
useful in various areas of physics, ranging from the classical description of the Higgs sector 
of the standard model (see for a review) to recent ideas in string theory (see |^ and 
references therein). 

This last example involves what we will call NonCommutative Yang-Mills (NCYM) the- 
ories in the sequel and can be thought of as a generalization of non-abelian gauge theories, 
whose gauge symmetry and interactions involves the noncommutative nature of the coor- 
dinates. A first example of such a theory appeared almost ten years ago, when Connes 
and Rieffel developed classical two dimensional Yang-Mills theory on the noncommutative 
torus 1^. This idea has also been generalized to higher dimensions 
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This naturally raises the question of the quantization of such theories, which has been 
tackled, at the one loop order, on the tori in [0 and on noncommutative M'^ in |Q and 0. 
Although these theories turned out to be non local, i.e. their interacting vertices involve 
trigonometric functions of the incoming momenta, it turns that the one-loop behavior is 
quite similar to the standard non-abelian case. This relies on an older work of Filk 
who proved that the trigonometric factor of any planar graph (in a sense to be defined 
below) does not involve trigonometric functions of the internal momenta flowing into the 
loops, thus exhibiting the same divergence as the standard theory. 

However, this is not true for non planar diagrams whose trigonometric factor does 
involve a phase depending on the internal momenta. Obviously such a phase softens the 
ultraviolet behavior of the corresponding diagram and it has been conjectured by several 
authors that such a diagram is in fact finite |]^, |TT], 



Nevertheless, it has been pointed out that this is not always the case |jT^. Indeed, if the 
non planar diagram contains some special kind of non planar subdiagrams whose standard 
degree of divergence is strictly greater that zero, which is the case in a scalar field theory, 
the small momentum behavior of these diagrams yields a new kind of infrared divergence 
intimately tied up with the non locality. 

This short paper is devoted to a survey of this problem in the simple case of multiloop 
corrections to the ghost propagator involving only nested and disjoint subdivergent one 
loop corrections to the ghost propagator. In the following section we shall briefly review 
the problem raised in |]13[ and then we shall present an explicit computation of the corre- 



sponding diagrams, postponing a complete probe into the renormalization of NCYM theory 
to a future publication. 

2 Small momentum singularities induced by non planar diagrams 

Before entering into the details of NCYM theory, let us recall that the noncommutative M'^ 
is the algebra generated by D hermitean elements with commutator = —2i6^jj, 

where O^y denotes a real antisymmetric matrix which we will assume to be of maximal rank 
for convenience. Furthermore, one introduces Fourier modes U{k) = e^'^'^, with k-x = k^x^. 
We will always think of a smooth and at infinity rapidly decreasing function as a Fourier 
transform 

/= fd'>kf{k)U{k) , 



where k f{k) is itself a smooth and rapidly decreasing function on standard M^. The 
commutation relations of the coordinates endow the algebra with the star product 



f^e9:= d''kd^lf{k)g{l)U{k)U{l) , (1) 



which yields 



with 6{k, I) = Ofj^^k'^l". Finally, this algebra is equipped with the analogue of an integral 
defined as 



// := /(O) (2) 
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and partial derivatives 



dj:= Jd^kik,f{k)U{k) , (3) 

which satisfy most of the properties of their commutative counterparts: positivity and def- 
initeness of the integral, Leibniz rule, commutativity of partial derivatives and integration 
by part, together with the tracial property of the integral 

f-^eg = jg-^ef, 

which proves to be fundamental in the construction of gauge invariant theories. 

At that point, two additional remarks are in order. First of all, let us notice that 
the matrix 6^1, explicitly breaks Lorentz invariance (or its euclidian counterpart), which 
is reduced to the transformations commuting with 6, whereas translational invariance is 
preserved. Furthermore, as 6^^^, is dimensionful, it also breaks scale invariance already at 
the classical level, for instance, in the case of four- dimensional NCYM theory or for a 
two-dimensional scalar field theory. 

From now on, one easily constructs scalar field theories, like 0^, whose euclidian action 

is 

f /I TTi^ g \ 

S[(f)]= / i^-d^4>'kgd^4)+ —4)'ke^+ —^-ke4>'^9^-^e^) , (4) 

or the NCYM action 

S[A^] = -^jF,,^eF''\ with (5) 
F^^ = df,A^ - d^A^ + g{Af, -kg A^ - A^ -kg A^) . 

Because of the tracial properties of the integration, the latter enjoys invariance under 
noncommutative gauge transformations 

SxAf, = g{X -kg A^ - Af, -kg A) - df,X . 

Perturbative quantization of these theories is easily performed within a formal functional 
integral point of view: The quadratic parts of the actions are equal to their commutative 
counterparts, whereas the interactions are non local and exhibit trigonometric functions 
of the incoming momenta of the interaction vertices. Thus, the total contribution of any 
Feynman diagram can be written as the product of a rational function by a trigonometric 
function. Because trigonometric functions are bounded, the standard rules of powercount- 
ing are unchanged so that Weinberg's convergence theorem remains valid. 

It has been shown that for planar diagrams the trigonometric function is indepen- 



dent of the internal momenta so that the Feynman integral reduces to the one encountered 
in a commutative field theory. For non planar diagrams the situation is more involved and 
we mainly have to distinguish two cases: whether the non planarity results from crossing 
of internal lines (i.e. from the non planar character of the amputed diagram), which we call 
type I diagrams, or whether it solely comes from crossing of internal lines whith external 
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lines (type E). The latter case is much more tricky because the phase vanishes when the 
corresponding external momenta satisfy some particular relation. 

The corresponding Feynman integral has been evaluated in [O within Schwinger's 



regularization scheme. It turns out that a type I non planar diagram whose powercounting 
subdivergent non planar diagrams are all of type I will be convergent, the corresponding 
singularity when the Schwinger parameters goes to zero being removed. If the diagram is 
of type E, but does not contain any subdivergence of type E, it is non singular, except for 
some exceptional values of its external momenta. 

Most of the trouble comes from the insertion of type E non planar sub divergences. 
Indeed, the latter correspond to Feynman integrals of the type 

j (T^'kRih, kn,Pi, . . . ,Pn) e'(^('=^'^i)+-+''('=-^")), (6) 

where fci, . . . , fc„ are the independent internal momenta, pi, . . . ,Pn are the external mo- 
menta and Pi, . . . ,Pn are linear combinations of the external momenta. The rational func- 
tion R{ki, . . . , kn,pi, . . . ,Pn) is responsible for the subdivergence. 

When all the Pj's which couple to internal momenta belonging to divergent integrals 
do not vanish, the corresponding integral can be considered as finite, being regularized by 
the oscillatory factor. Indeed, a computation with a cut-off introduced within Schwinger's 
parametric formula yields such a finiteness. However, whenever one of the Pj's coupled to 
a divergent loop integral vanishes, then the corresponding Feynman integral is just a usual 
divergent loop integral and yields a singularity. When inserted into a larger diagram, this 
could create some trouble when integrating over momenta approaching the subspace Pi = 0. 
The question is how fast the divergence appears compared with the smoothening property 
of the integration measure. A quadratic divergence seems to destroy the renormalizability 
13| whereas a logarithmic divergence could be harmless. This conjecture is supported by 



our simple example. 

In a mathematically more satisfying manner, one can also consider such an integral 
as a well defined distribution which is nothing but the Fourier transform of the rational 
fraction R. However, we still have to face a problem when inserting the distribution into 
a larger Feynman diagram. In particular, the example described in |jl3| corresponds to a 



distribution which is nothing but the Feynman propagator and the infrared troubles when 
p ^ are quite similar to the usual small x singularities encountered in QFT. 

Finally, let us point out that we did not encounter such a problem when quantizing 
NCYM theory on a torus . The construction of the latter is similar to that of NCYM 
on except for the quantization of the momenta appearing in the Fourier transform. As 
a consequence, there is no singularity when p — >■ since such a limit cannot be taken. 
However, we noticed that a one loop non planar diagram with external momenta p exhibits 
an extra UV singularity of the type 6{p)/e. Quite miraculously, all these singularities turn 
out to cancel, leaving us with finite one loop renormalized correlators. 

In the next section, we shall evaluate the non planar contribution to some of the one 
loop corrections to the ghost self energy using Bessel functions, showing that they do not 
lead to any singularities when inserted into larger diagrams. We shall use the Feynman 
rules for NCYM theory without deriving them, the latter being obtained by replacing the 
structure constant of non-abelian gauge theory fabc by 2isin6'(p, g) M. 
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3 A simple example. One-loop calculation 

Our goal is to compute the following 1-loop correction to the ghost propagator in NCYM 
theory: 

k 

p ) k+p \ p (7) 



Wavy lines represent gluons and straight lines ghosts. The Feynman rules derived in 
and (for the noncommutative torus) |]^ lead to the integral 

/i = |rf'*4g^R(-p-*)''^(^;..-(l-a)^)(|;j^(-p)"'sta9(*.p) sme(-k,k+p) 

We work in a = 4 dimensional euclidian momentum space with metric g^i, = 6^i, and 
use obvious abbreviations such as pk = g^^p^k'' . Using Feynman parameters 

V{r+s) x'-Hl-xy-^dx 

(9) 



A'-B' T{r)T{s) Jo {Ax + B{l-x)Y+' 
we obtain 

' (p'^+apk) 2{l-a){l-x) {pky 



I-^ = -Ag'^n jd^k sin^ e{k,p) j^dx(^ 



{k'^+2pkx+p'^xy {k'^+2pkx+p^x) 



In the denominator we write qk instead of pk so that we reproduce the /c's in the numerator 
by differentiation with respect to q: 



h = -^g^h Jd^k sm^9{k,p)(^p^ 



dx ,, d f adx 
pM. 



(k"^ + 2qkx + p'^xy dq^^ Jq 2x{k'^ + 2qkx + p'^x) 
d d {l—a){l~x)dx 



p^p"- 



dqfj- dq^ Jq Ax'^^k'^ + 2qkx + p'^x) 
Using = dtt"'~^ e~^^ we rewrite the integral into 

. 2. f\ f 2 01 B (l-a)(l-x) „ , d d \ ^^, 

K[t,p,q,x] := j rf^A;e-*('='+2«'=^+f'")sin2^(A;,p) . 
Developing the sine into a Fourier series we obtain for the kernel 

K[t p q X] = J d^k (^^Q-t('<^+1^)^-t^(p'^^-l'^^^) _ lQ-t(k+qx+ie{p)/t)^-t{p'2x-q^x^)-2ixe{p,q)-pop/t 



— lp-t(k+qx-ie{p)/ty -t{p'^x-q-^x^)+2ixe{p,q)-pop/t\ 
4^ J 
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where 9{p)'^ := O^'^Pa and pop := g^yO^'^Pa 6"^pp = {6{p)Y. Then it is easy to perform 
the Gaussian integration: 

2 

K[t, p, q,x\ = ^ (^^-t{p^x-q^x^) _ ^~t(p^x-q^x^)-pop/t 2x0(p, q)) . (10) 
We can now perform the differentiations with the following result: 

I-^ = Tx"^ g'^n I dx I dt{^p^{?,ax-a-l-x) + 2{p^f{l-a)x'^{l-x)) x 

X (^Q-^P^i'^-^) — Q-tp^xil-x)-pop/f^ ^^^^ 

g^n dx dt{-{3ax-a-l-x) + 2{l-a)x) (e^* - e-t-x{i-x)p^pop/t^^ _ 
Jo Jo t 



2 2 2i 
p IT " ^ 



The integration over t diverges logarithmically. We define the projection tp onto the diver- 
gent part by 

tl{I,) := p\^g^h ['dx r -{3ax-a-l-x)e-' = -p\^g^h^ H "—^ . (12) 
Jo Jo ^ 2 Jo ^ 

The projection to the convergent part is 

R,:={l-tl){h) ^ ^ (13) 

= p^7rV^((l-«) +2(l-a)x)e-^-"(^-")^'p°^/*). 

Instead of introducing a cutoff as in [|l^ we prefer such a momentum subtraction before 
performing the divergent integral, in analogy to the BPHZ scheme, as the notation 
for the projection indicates. We believe this is advantageous for the renormalizability 
proof to all orders based on Zimmermann's forest formula. Please notice the difference 
between the power counting degree 2 entering the forest formula and the actually only 
logarithmic divergence in (|T^. This is crucial for the insertion as sub divergences and 
gives the reason why we will obtain local counterterms to all orders whereas there are true 
quadratic divergences in the scalar theory in |]13|. It should be not difficult to extend such 
a subtraction scheme to the entire NCYM theory. It is however important to define 
as the projector onto the strictly divergent part of an integral in order to produce local 
counterterms. This assumes one can prove that all divergent integrations give rise to such 
local terms, as we will do in this paper for repeated insertions of the ghost propagator. 

The x-integration in (|TB|) yields Jq dx ^ e"*""^^^"^^/* = ^ dx e~*~"^(i~^)/* for any a so 
that 

i?, =p2^V;i((i-«)+ fdx nt(^-(i-«))e-*-^(i-^)p'^°p/*). 

^ Jo Jo 2t / 

The t-integration leads to Bessel functions: 

Ri = p^Tx^g^h{l—a){l — j dx 2^/x{l—x)p'^p o p Ki[2^/ x{l—x)p'^p o p]j 

+ p^TT V/i (3-a) / dxKo[2^/x{l-x)p'^pop] . (14) 
Jo 
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The Bessel function Ko[y] diverges logarithmically to +00 for y — > and converges expo- 
nentially to for y — > +00. Hence, for any exponent r > there exists a number c° > 
such that 



Ko[y] < c'Jy^ 



V < y < 00 . 



(15) 



This will be proven algebraically in the Appendix. Graphically the situation is sketched in 
Figure 1. The function yKi[y] approaches 1 for y — > and converges exponentially to 



3 
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2 

1.5 

1 
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Figure 1: The Bessel function Koly] (lower 
graph) is for any < y < cxo bounded 
by the function c°/y'' (upper graph). The 
situation is shown for r = 1 and = |. 
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for y — > +00. It is nevertheless convenient to regard it as KQ[y] before: For any exponent 
r > there exists a number > such that 

yK^[y]<cl/y' V < y < 00 . (16) 

The corresponding graphic is shown in Figure 2. Now the x-integration is easy to perform. 



Figure 2: The function yKi[y] (lower 
graph) is for any < y < 00 bounded 
by the function cj/y*" (upper graph). The 
situation is shown for r = 1 and c\ 
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It converges for < r < 1, which means that infrared divergences are absent (for non 
exceptional momenta). 

We restrict ourselves to the case where the rank of the tensor g^vd^^^O^^ equals the 
space-time dimension (maximal noncommutativity) . Then there exists some parameter 
mp of dimension of a mass (the 'Planck mass') such that 



po p > 



p 



ml 



p 



ml 



This yields the estimation 

R,^p'nVh{{l-a) + 0[pXa){'^ 



(17) 



(18) 



where -P"(a) is a polynomial of homogeneous degree n in (1— a) and {3— a) with coefficients 
of order 1. 
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4 Higher loop order calculation 



Now we insert n of these 1-loop propagator corrections into a propagator correction, giving 
an n+1 loop diagram: 



n+1 




(19) 



Using d^) and the Feynman rule for the ghost propagator, the corresponding integral is 
of order 



+ apk 
¥{kTW 



(l-a) 



{pkf 



r) 



X 



3=0 



—a 



k^k^ik + py^ 

^ (20) 



(k + p) 



with -P°(a;) = 1. Introduction of Feynman and Schwinger parameters as before and use of 
xT{x) = r(a;+l) leads to 



= _4/^(_vrV/i)"^(l+jr)P/(a) (l-a)"-^m^^" d^k sin^ e{k,p) x 
j=o 

x^'^ijp'+apk) {2+ir){l—a){l—x)x^'^{j)kY 



Jq \{k'^+2pkx+p^xf+^'' 
^ r(l+jr) 



(A;^+2pA;x+p^x)^+-^'' 



^m^^'^ X 



rfx dt^-x^'t^+^'p^ + |ax-'"-^t^>'^ 



+ i(i-«)(i-^)x^'^-Vw^^)i^[t,j>,g,x] 



2 2 
-TT ^ 



X j^dx dt(^{3ax-x-a-l)xn''''y + 2{l-a){l-x)x^'+h^%py^ x 



1 

r(l+jr) 



■^rrip X 



-TC' 



r(l+jr) V / Jo (l-a;)-"^ 



X / dt{{3ax-x-a-iy'-' + 2{l-a)xf'') (e"* - e-*-(i-^')pVp/t) _ 



(21) 





The only divergent integral is for j = the projection 

>1 POO 



tl{In+i) := -(-7rV^)"+V(l-a)" / ^^a; / -(3«x-x-a-l)e- 



(-7r2/;.)"+y(l-«)' 



3-a f^dt 



— e 



(22) 
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The convergent part can be evaluated to 

Rn+l = (1 - ip)(-^n+l) 

= -p2(_7rV^)"^^(l-a)''^^(l - ^ dx 2 v^x(l-x)p2p o p Ki[2^x{l-x)p'^p o p] 

-j92(_vrV^)""^^(l-a)"(3-a) [ dx Ko[2^/x{l-x)p^p o p] 

Jo 

+ A~^Vnr^ ± (l-a)"-(3-a) i^fy (23) 

" 1 2 ■ 

p1 j poo 

The last two lines range from zero (for p = oo) to minus the value of the third last line 
for p = 0. Thus, in our estimation we have to neglect the last two lines. The remaining 
integrals over Kq and Ki are familiar to us, see Figures 1 and 2, and we choose the essential 
exponents in ([15|) and (|l^) to be (n+l)r instead of r. Then we arrive at 



n+l 2 




But this was precisely our starting point we inserted into (||) to obtain (pO]). Hence, 
(P^) provides the structure of any renormalized n+l loop graph made of ghost propagator 
corrections. The counterterm of such an n-loop graph is given by (^), and we see explicitly 
that the Feynman graphs made of nested 1-loop ghost propagator corrections 



(25) 



are renormalized by local counterterms for any order in %. Here, locality means that the 
momentum dependence of the counterterm and the kinetic part of the ghost action are 
identical. In order to renormalize an n-loop graph (to avoid infrared divergences) the 
critical exponent has to be chosen < r < 1/n. 

The essential step in this proof was the observation that the noncommutative Feynman 
graphs under consideration evaluate to Bessel functions, which can be estimated by a power 
law. It seems plausible that any Feynman graph of noncommutative Yang-Mills theory 
evaluates to Bessel functions, and applying the same techniques it should be possible to 
show that local counterterms suffice to renormalize this model. 

Appendix: Proof of Eq. (|T5|) 

We prove that for each r > there is a number > such that 

Coix) ■= — > Ko(x) yO<x<oo. 
x'^ 
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The Bessel function Kq{x) is one of the two solutions of the differential equation 



xKq + K'q — xKq = , < X < oo . 
It is however more convenient to consider the function 



(26) 



K{x) := ^/^Ko( 



x] 



and compare it with 



C{x) := V^Coix) 



1 — Ax"^ 



1 - 4r^ 
C" + — 4^C = 



The derivative of the Wronskian W{K, C) := K'C - C'K is 



W = ^—^KC 

x^ 



W > Q for X > r 
W <Q for X < r 



The asymptotic development shows that for the solution Kq of (|26|) one has W{x) < for 
X —>■ oo and W{x) > for x — > 0. Therefore, there is only one zero of the Wronskian, at 
illustrated in Figure 3. 



0.3 
0.2 

0.1- 

-0.1- 




Figure 3: The Wronskian W{K, C) has an 
extremum for x = r and a zero for x = x^. 
The situation is shown for r = 1 and Cr = 1. 



We choose the normalization 
so that W{xr) = and K{xr) = C{xr). Then we can integrate 



r-l 



^ K'(x) , '•^ ^' 

ax < 



Xr 



K{x, 
K\x) 
K{x) 



C'jx) 

Xr C{X) 



dx 



dx > I —r^—dx 



X C{x) 



K{x) < C{x) 
K{x) < C{x) 



X > Xr ^ W{x) < : 
X < Xr ^ W{x) > : 
This finishes the proof. 
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